Quantum Mechanical Search and Harmonic Perturbation 
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Perturbation theory in quantum mechanics studies how quantum systems interact with their 
environmental perturbations. Harmonic perturbation is a rare special case of time-dependent per- 
turbations in which exact analysis exists. Some important technology advances, such as masers, 
lasers, nuclear magnetic resonance, etc., originated from it. Here we add quantum computation to 
this list with a theoretical demonstration. Based on harmonic perturbation, a quantum mechanical 
algorithm is devised to search the ground state of a given Hamiltonian. The intrinsic complexity 
of the algorithm is continuous and parametric in both time T and energy E. More precisely, the 
probability of locating a search target of a Hamiltonian in iV-dimensional vector space is shown 
to be 1/(1 + cNE~ 2 T- 2 ) for some constant c. This result is optimal. As harmonic perturbation 
provides a different computation mechanism, the algorithm may suggest new directions in realizing 
quantum computers. 
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INTRODUCTION 



Quantum physics can in principle speed up solving the 
unsorted-database search problem with a quadratic im- 
provement over classical algorithms, as was first demon- 
strated by Grover [ij. This problem was originally for- 
mulated as to identify a target item in the fewest queries 
to a black-box database. An important reformulation 
by Farhi et al. 0, Q phrased the problem as to search 
the target state with some special eigenvalue of a given 
Hamiltonian TL, which represents the database. In partic- 
ular\ constant and adiabatic perturbations were proposed 
in 0] and respectively, for quantum search. These 
methods yield the same quadratic speed-up as Grover's 
construct : S\. Moreover, the adiabatic computation is 
equivalent to standard quantum computation (in terms 
of unitary transformations) (f|. Although some physi- 
cal implementations have been demonstrated in realizing 
quantum search algorithms, e.g. 0, HQ, they may not 
be scalable to solve large problem instances without some 
fundamental breakthroughs [lo| . Searching alternative 
computation models may suggest new ways of building 
quantum computers. 

Despite the success of the constant and adiabatic per- 
turbations in quantum search, the applicability of pertur- 
bations based on fast time- varying Hamiltonians remains 
an open problem. This paper exploits harmonic pertur- 
bation for quantum computation. A new computation 
model is proposed, inspired by the well-studied harmonic 
perturbation of two-state systems in quantum mechan- 
ics. By preparing a system in one or the other of its two 
states initially, such perturbation induces an absorption- 
emission cycle (the phenomenon of periodic oscillation 
of the probability for the system being found at one of 
its two states) at the resonance condition. More specifi- 



cally, consider a two-state physical system of Hamiltonian 
H = £u|l)(l| + E 2 \2)(2\, with E 2 > Ei, in a sinusoidal 
potential V(i) = 7e Mt |l)(2| + T e~^*|2) (1|. The state 
evolution of the system is governed by the Schrodinger 
equation 



'' dt 



(1) 



with ip(t) = ci(t)\l) + c 2 {t)\2) and 7i(t) = H + V(t). 
For the initial condition ci(0) = and c 2 (0) — 1, the 
respective probabilities of finding the system in states 
|1) and 1 2) are of exact solutions (Rabi's formula; see, 
e.g., El) 



|ci(*)| 2 = 
|c 2 (£)| 2 - 
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7 2 + (w - W2i) 2 /4 



sin 2 (fit) 



(u-U 2 l) 2 



(2) 
(3) 



At 



where ui 2 i — {E 2 — Ei) and ft = y 
resonance, to = u> 2 i, the probability of finding the sys- 
tem in the ground state |1) oscillates with period 71-/7, 
and reaches 1 at time (2k + l)n/(2j), k = 0,1,2,.... 
This phenomenon reveals the potential usefulness of har- 
monic perturbation in searching the ground state of a 
given Hamiltonian even with arbitrarily multiple states. 

Intuitively harmonic perturbation to a general 
multiple-state system may be used as a mechanism 
searching the ground state of the corresponding Hamilto- 
nian. Knowing the energy gap between the initial state 
and ground state, one can apply a perturbation at the 
resonance frequency to induce a probability oscillation 
such that the system swings between these two states. If 
a measurement is performed at the right time, the system 
situates definitely in its ground state. Based on this intu- 
ition, we devise a quantum search algorithm. An analysis 
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shows that, given a Hamiltonian in TV-dimensional vector 
space with energy gap E between its ground state and the 
other TV— 1 excited states, our algorithm finds the ground 
state in time T with probability Pr = l/(l + cNE- 2 T- 2 ) 
for some constant c. Therefore, to search the right tar- 
get with a constant high probability independent of N, 
the time-energy product complexity matches prior known 
O(viV) result [12j of other quantum search algorithms 
[E S Sj- Note that, due to the time-energy duality, it 
is not meaningful to speak about only time or energy 
complexity, regardless of the other. 



Algorithm 1: Naive Ground-State Search 
input: a Hamiltonian TL 
output: the ground state of TL 
begin 

01 measure state, yielding say \gi). 

02 if E gi = return \ gi ) 

03 apply harmonic perturbation V[ ffll7l B] (t). 

04 measure state at time 71/(27), yielding say 

05 return |<?2) 
end 



52}. 



FIG. 1: The procedure of a naive quantum search. 



QUANTUM SEARCH ALGORITHMS USING 
HARMONIC PERTURBATION 



We formulate the database search problem as follows. 
Given a Hamiltonian TL = 'EjL 1 Ej\j){j\, we are asked 
to find state \g) such that E g is the minimum among 
Ej's, i.e., \g) is the ground state. In the sequel, we as- 
sume any state \j) is a configuration of n binary digits; 
thus, N = 2 n . Moreover, we assume that the state of 
the underlying physical system is measurable such that 
its n-bit configuration (e.g., spin orientations of spin- 1/2 
particles measured along some axis) is completely deter- 
mined, and that its corresponding energy under TL can 
be observed thereafter. Unless otherwise stated, we shall 
focus on the energy distribution of Grover's search prob- 
lem, and assume one out of the N states is the ground 
state of energy and the other N — 1 states are excited 
states of energy E > 0. In our algorithm, we consider E 
is adjustable in analyzing energy complexity. 

Below we show the new application of harmonic per- 
turbation in quantum database search. In principle, if 
the perturbation potential is designed properly, knowing 
the initial (excited) state and energy gap of a given N- 
state system, one can apply the resonance frequency to 
induce an oscillatory transition almost solely between the 
initial state and the ground state, similar to the two-state 
case. Measuring the system at the right time achieves the 
highest probability of locating the ground state. How this 
peak probability is related to N is our main concern. To 
gain an insight on the optimality limit of quantum search 
using harmonic perturbation, we begin with a simple trial 
(in Figure [T]) and then proceed with an optimized proce- 
dure (in Figure |4]). 



A trial algorithm 

Figure[T]sketches a (non-optimal) quantum search pro- 
cedure. The algorithm starts with a measurement in 
the state basis to enforce the underlying physical sys- 
tem collapsing to some state, say \gx). If the correspond- 
ing eigenenergy of the input Hamiltonian TL in state |<7i) 
equals 0, the algorithm has found the target and returns 



\gi) immediately at Step 2. Otherwise, harmonic pertur- 
bation is applied using the sinusoidal potential Vu )7)(1 ,] (t) 
with 



(p|v b - 7>tl> ] (*)|?> = 



if q = j and q 7^ p 
if p = j and p =/= q 
otherwise 



for indices p, q = 1, . . . , N. By replacing index j with g\ 
and letting u> — E, the perturbation V\ g ^ tl m(t) at Step 3 
induces an oscillatory probability for the system swing- 
ing mainly between state \g\) and the unknown ground 
state. The algorithm measures, at Step 4, the state of 
the system at time ir/(2y), when the system situates in 
the ground state with the highest probability. 

We analyze the condition under which this peak prob- 
ability is independent of the effect of N and is close 
to 1. For N = 2, the returned state \g2) of the al- 
gorithm is the ground state with certain. However, it 
is not the case for N > 2. To see why, we solve the 
Schrodinger equation of the iV-state system with Hamil- 
tonian TL(t) — Tl+V[g in .E] (*)• To simplify the discussion, 
assume without loss of generality that |1) is the ground 
state and \gi) = |2) is the initial state. Let c/-(i) denote 
the probability amplitude of state \k) at time t. Then the 
original N first-order differential equations from Equa- 
tion |T]) can be reduced to three due to the equivalence of 
C3(i), Ci(t), • • • , cjv(i)- Hence, at resonance u> = E = ojr, 
the reduced equations in terms of bk(t) = e lEkt Ck(t) are 



ib\(t) = 762(4) 

ib 2 (t) = ybi(t) + (N-2)ye- iu «%(t) 
ib 3 (t) = 7e^'62(t) 



(4) 
(5) 
(6) 



To solve 1 61 (t) 1 2 and thus |ci(i)| 2 , apply Laplace trans- 
form C on these equations and solve for -Bi(s) = C{b\ (t)} 
with initial conditions 62(0) = 1 and 6fc(0) = for k 7^ 2. 
We derive 



Bx(s) 
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1 



( s2 +7 2 ) (l + j^AO 



(7) 



where A 1 = (TV — 2)s/ (s + iujR). From the inverse Laplace 
transform of -Bi(s), the exact solution of Ci(i) can be 
derived. For N — 2 and thus Ai = 0, ci(i) reduces 
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FIG. 2: For 7 = 0.01, lor = 10, the figure plots |fei(t)| 2 of Al- 
gorithm 1 under various N and t. It shows that max* \bi(t)\ 2 
decreases and |&i(t)| 2 oscillates faster as N gets larger. 



FIG. 3: For 7 = 0.01, the figure plots max t |6i(i)| 2 of Al- 
gorithm 1 under various N and lor. It shows that letting 
lor oc N maintains max* |&i(t)| 2 at some constant value. 



to Equation ([2]). For general N > 2, we omit listing 
the sophisticated expression of exit) as the trial pro- 
cedure is not our final destination. Nevertheless, the 
numerical simulations as shown in Figures [2] and [3] re- 
veal the performance of Algorithm 1 with respect to 
N. An analysis suggests that, for a fixed constant 
7 (among other possibilities), to maintain a constant 
peak probability, Pr = max t |ci(i)| 2 , we need ojr oc N . 
Thus, the algorithm has time complexity 0(1) due to 
the fixed perturbation amplitude 7, and has energy com- 
plexity 0(\\H\\ 2 + ||V|| 2 ) = 0{uj r + jVN) = 0{N) for 
lor G O(N). This constant complexity in time and linear 
complexity in energy can be explained in the s-domain 
by observing that N and lor in Ai of Equation are of 
a first-order relation. That is, by maintaining N/lur = e 
for some small constant e, the effect of Ai is negligible 
and |ci(t)| 2 s» sin 2 (jt). Hence, the algorithm achieves 
the same linear resource complexity as the classical algo- 
rithm for database search. 

As a digression, we noticed that the recent work [l3| 
presented an algorithm on quantum search with reso- 
nance, where the perturbation is very similar to Vu )7 (t) 
of Algorithm 1 and thus the algorithm is unlikely to 
be superior to the classical counterpart. The miscon- 
ceived quantum improvement may be due to the igno- 
rance of the time-energy product complexity as well as 
due to over-simplified analysis [lij]. Moreover, the formu- 
lation assumed that the eigenenergy of the search target 
is unique and known a priori. However, under this as- 
sumption, the energy eigenstates could be pre-computed; 
the target eigenstate can be known once the eigenenergy 
is specified without even resort to quantum search. (In 
contrast, we assume that the excited energy eigenstates 



Algorithm 2: Optimized Ground- State Search 
input: a Hamiltonian TC 
output: the ground state of TL 
begin 

01 measure state, yielding say \gi). 

02 if E gi = return \ gi ) 

03 apply harmonic perturbation V° gi 7 e}{^)- 

04 measure state at time 77/(27), yielding say \g2). 

05 if E g2 — return \g-2) 

06 apply harmonic perturbation V^ g2 „ B j (t) . 

07 measure state at time 77/(27), yielding say |<?3). 

08 return ^3) 
end 

FIG. 4: The procedure of an optimized quantum search. 



are degenerate and there is another observable commute 
with the Hamiltonian that further determines the states.) 



A refined algorithm 

Based on the insight from the s-domain analysis, we 
obtain an optimized algorithm. Figure 2] sketches a re- 
fined quantum search algorithm similar to that of Fig- 
ure [T] It differs from the procedure of Figure [T] mainly in 
an additional iteration, and in the applied perturbations 
0) and v y, 7)W ](*)> where 



VP. 



(P|V&, 7JU ](*)|9> 



ye 



7e" 



J * if (p — j q and q even) or 

(q = j 7^ p and p odd) 
-tut jf (q = j 7^ p and p even) or 

(p = 3 7^ <Z an d q odd) 
otherwise 
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for indices p, q = 1, . . . , JV and Vg.^t) = V^_ u] (t) 

( Le -' V [°, 7 ^]M and V [i : 7^](*) are conjugate to each 
other). For instance, V° 2 w i (t) in a 4 x 4 matrix reads 

/ 7e Mt \ 

7e Mt 
V 7e- Mt / 

Because Vg w i (t) and V£ w i (i) induce noticeable prob- 
ability oscillation only when the ground state, say \g), 
situates at an odd (i.e. g odd) and even (i.e. g 
even) position, respectively, Algorithm 2 requires one 
more perturbation-and-measurement iteration than Al- 
gorithm 1. If the ground state, not returned in Step 2 
of Figure [H situates at an odd (respectively even) posi- 
tion, it will be returned at Step 5 (respectively Step 8) 
almost for sure. The algorithm can be repeated to fur- 
ther enhance its correctness probability. (Notice that, in 
preparing V° and V e , we need to specify the initial state 
for parameter j. Accordingly the initialization measure- 
ments at Steps 1 and 4 determine the parameters j of 

v [°-,7^]( < ) and v bV7>">](*)' res P ectivel y- That is wh y the 

perturbation V e of Step 6 depends on the measurement 
result 1 52) of Step 4. Alternatively one may discharge 
this dependence by fixing j to some specific state.) 

We analyze the condition under which the probability 
for Algorithm 2 finding the search target is independent 
of the effect of N and is close to 1 . To simplify our discus- 
sion, again we assume without loss of generality that |1) is 
the ground state of TL and the first measurement in Step 1 
of Figure |4] yields \g\) — |2). To compute the probability 
that the ground state |1) is correctly returned in Step 5, 
we solve the Schrodinger equation of an iV-state system 
with Hamiltonian TL(t) = TL + V° 2 7 E j(t). Let Ck(t) be 
the probability amplitude of state at time t. Then 
the original N first-order differential equations can be 
reduced to four due to the equivalence of odd coefficients 
cs(t), . . . , cjv-i(t) and the equivalence of even coefficients 
C4,(t), . . . , C]y(t). Hence, at resonance lu = E = lor, the 
reduced equations in terms of bk(t) = e lEkt Ck(t) are 

ib\{t) = 762 (*) (8) 
ib 2 (t) = 761 (t) + 

{N 2 2) l(e- lURt b 3 (t) + e iu *%(t)) (9) 

ib 3 (t) = 7e^« t 6 2 (t) (10) 
ib 4 (t) = ie - iURt b 2 (t) (11) 

To solve 1 61 (t) 1 2 and thus |ci(i)| 2 , we apply Laplace trans- 
form C on these equations and solve for -Bi(s) = C{b\{t)} 
with initial conditions 62(0) = 1 and fcfc(0) = 0, k 7^ 2. 
We derive 



where A 2 = (N— 2)s 2 /(s 2 +luji). Taking inverse Laplace 
transform £~ 1 {Bi(s)} and assuming N ^> 1, we get 

h( t ) « sin ( m t ) (13) 

By Equation (| 13[> . the peak probability Pr of finding the 
search target equals max t |ci(i)| 2 = 1/(1 + N^ 2 ujji^ 2 ); 
the period of the probability oscillation of |ci(t)| 2 is 
t = TTy/l + A r 7 2 cjfl _2 /7. Since on average the search 
target can be found by running the algorithm 1/ Pr times 
each of which takes time r, the total time complexity is of 
r/ Pr S 0((1 + N^lur' 2 )^ /j). On the other hand, the 
energy complexity is of 0(||i?||2 + || VH2) = 0(ujh+jVn). 
It can be verified that maintaining a constant Pr achieves 
the tightest upper bound of resource complexity. As a 
result, when Pr is maintained as a constant (as we shall 
assume in the sequel), 7 is the only parameter affecting 
time complexity T. Letting E denote the energy com- 
plexity, we can write Pr = 1/(1 + cNE~ 2 T~ 2 ) for some 
constant c. In essence, the time-energy product complex- 
ity is of 0(VN). 

The foregoing analysis assumes that the ground state 
is in an odd position. Suppose that the ground state is 
rather in some even position. Then V° gi 7 E At) of Step 3 
does not induce a state oscillation (and thus \g 2 ) is the 
same as \g\) with high probability) while Vf g2 of 
Step 6 does. In this similar analysis holds to de- 

rive the probability that the algorithm returns a correct 
answer at Step 8. 

Since Algorithm 2 achieves the 0{\/~N) time-energy 
product complexity, it is optimal according to the anal- 
ysis of [2|. Below we provide yet another proof using 
s-domain analysis. It may provide an intuition why Al- 
gorithm 2 outperforms Algorithm 1. Also we hope that 
the technique of s-domain analysis may become a pow- 
erful tool for the complexity analysis of other quantum 
algorithms. As least in our cases, s-domain analysis is 
more intuitive than time-domain analysis. 

THE OPTIMALITY LIMIT 

Comparing Algorithms 1 and 2, one might think that 
the performed modification might seem suspiciously irrel- 
evant as it only affects states |3), . . . , \N) in our example. 
However, this is not true because it in fact affects |2) by 
Equation © and thus |1) by Equation ©. (Note that 
Algorithm 2 on perturbation V°, and also V e , applies 
both e luJt and e~ luJt simultaneously, rather than sequen- 
tially. It differs from applying e tuJt on some states for 
one search and then applying e~ luJt on other states for 
another search. Otherwise, Algorithm 2 would be equiv- 
alent to performing Algorithm 1 for several times and 
has no improvements.) Although the effect may seem 
obscure from the time-domain Schrodinger equations, it 
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is apparent when solving them in the s-domain. The 
two different phases in V° and V e of Algorithm 2 may 
interfere and simplify i?i(s). 

To understand the improvement of Algorithm 2 over 
Algorithm 1, compare Equations ([7]) and (fT2")) . We see 
that the undesirable effect of large N in A2 is nullified by 
ujji while that in Ai is nullified by lor. As a consequence, 
for constant 7, Algorithms 1 and 2 are of energy complex- 
ities O(N) and 0(y/N), respectively. Ideally if we can 
construct a A such that the effect of N is cancelled out by 
lor for a larger exponent d, then Pr can be maintained 
as a constant with a lower resource complexity. Unfor- 
tunately, we will now show that d is at most 2, that is, 
no improvement is possible by introducing more differ- 
ent phases and/or amplitudes to the non-zero entries of 
V [°, 7 , W ]W and V [i, 7 , w ]W- ( Note that replacing the zero 
entries of w i (t) and Vt 7 w i (t) with non-zero elements 
introduces undesirable amplitude leakage to states other 
than the initial and target ground states. Hence we only 
need to consider modifying the non-zero entries.) 

Consider a general perturbation potential (a Hermi- 
tian matrix) with non-vanishing entries only in the row 
and column indexed by some initial state \j) similar to 
Vu 7iW i(i). These entries can have arbitrary amplitudes 
and frequencies. (However, one of the frequencies must 
equal lur such that resonance is possible.) Assuming 
without loss of generality |1) and |2) to be the ground 
and initial state, respectively, we solve the corresponding 
Schrodinger equations and can write 



B 1 (s) = 



1 



(* 2 +7 2 ) (i + I - 5 X_ 3 _A) 



with 7 being the amplitudes for entries |1)(2| and |2)(1|, 
and 



A = 



« 2 



IUJ 2 



(14) 



where index m is polynomial in n (thus m -C N = 2"), 
otj 's are positive real numbers, ojj 's are of the form ajU>R j 
for real constants aj and Xj , and u>j 7^ ojk for j 7^ k. (Note 
that dj is obtained from the product of some complex 
number and its complex conjugate, and thus is positive.) 



eralized problem that the search target appears more 
than once in the database (i.e-, there are multiple ground 
states in our setup). For [l|, the number of appear- 
ances needs to be known a priori for accurate search 



(unless more complicated quantum counting 14] is in- 
corporated); for 2], the method may also seem no easy 
extension. In contrast, adiabatic computation Q sup- 
ports multiple appearances of search target. However, it 
requires careful analysis to decide valid evolution speed; 
whether there is quantum speedup may depend on the 
instance to be solved. On the other hand, our approach, 
when applied to the generalized problem, may not al- 
ways yield quantum speedup. An analysis shows that 
a large number of ground odd-states (respectively even- 
states) for perturbation V° (respectively V e ) makes the 
multiplicities of 63 (i) and 64(f) in Equation ^ differ to 
some extent. Thus the advantage of Algorithm 2 with 
respect to Algorithm 1 disappears. It would be interest- 
ing to know if our approach can be extended to handle 
this generalized problem with definite quantum speedup. 
Also whether our method has unseen advantages over 
other methods remains further investigation. 

Even though our perturbation Hamiltonian looks more 
complicated than that of [2'], it does not directly imply 
our proposal is harder to implement. The implementa- 
tion issue is beyond the scope of our focus and poses 
a challenge to the experimentalists. Nevertheless, it at 
least suggests a new way of performing quantum search. 



CONCLUSIONS 

Motivated by developing a quantum search algorithm, 
we formulated a multi-state harmonic perturbation prob- 
lem and gave an exact analysis. This paper answered af- 
firmatively and constructively the open question whether 
fast time- varying Hamiltonians can be exploited in quan- 
tum search. The presented algorithm is optimal and 
achieves a quadratic speed-up over classical algorithms 
similar to prior methods [H, 0, 0] • Under this new com- 
putation model, we hope it may suggest different ap- 
proaches to the realization of quantum computers. 



Theorem 1 Let d equal the largest exponent in terms of 
ljji in the denominator of Equation J_?^[ ) minus that in 
the nominator. Then, d < 2 for all possible assignments 
to the constant parameters aj, aj, and Xj of A. 

Since d = 2 is achieved, the quantum search algorithm of 
Figure H] is optimal. 

DISCUSSIONS 

We compare various quantum search algorithms, 
namely, those of [H, 0, Q and ours, in terms of the gen- 
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Appendix 

Theorem Q] Proof: 

Expanding A of Equation (14)) , we have 



A = 



1 



-{s m + --- + 



(s + iui) ■ ■ ■ (s + iu m ) 



k^3 



We first show that, if d > 0, the largest exponent in 
terms of u>r in the denominator of Equation (|15|) must 



result from the product of all Wj's. Let a = x\ + ■ ■■ + 
x m ; then the product term ■ ■ 'W m equals auiif for 
some constant a ^ 0. For contradiction, suppose that the 
product term oj\u>2 ■ ■ ■ w m is not the only highest order 
term in u>r. Then there must exist at least one u)j — 
a,jUjR Xj with Xj < 0. Let a'ujR" for a 1 ^ be the product 
term of all lo^s with Xj > 0. Because a'ajR 17 must appear 
somewhere in both the denominator and nominator of 
Equation ([15]) . d must equal 0. Hence, to have d > 0, 
all ujj's must have Xj > 0. That is, the product term 
• ■ ■ u> m contributes to the largest exponent of ujr in 
the denominator. 

To have d > 2, the exponent of cor in the nominator 
of Equation ([15]) cannot be too large. Observe that, for 
the resonance condition to hold in harmonic perturba- 
tion, there must exist some Xk — 1. Moreover, there are 
more than one such x^s. Otherwise, the product term 
u>i • • • Wk-i^k+i ■ • • u> m cannot be cancelled out. Its exis- 
tence in the nominator of Equation [15] results in d < I, 
which violates the desired condition. Because of these 
constraints and Xj > for j = 1, . . . ,m, the coefficients 
of s and s 2 in the nominator of Equation (|15p must equal 
zero. That is, 



e^e n ^ 

must be satisfied. Multiplying (fl6|) by 
yields 



r 1 



£« 



n 



£« 




(16) 
(17) 

-- + £> 
= 0. (18) 



By recognizing that the second term of Equation {T 
equals by Equation (jTTj) , it is immediate that 



\ - rifc^j ^k 



(19) 



Again, multiplying Equation Q19[) by the product term 
luiu>2 ■ ■ ■ u>m yields 



3 k=£j 



(20) 



Since ay's are positive real numbers, the left-hand side of 
Equation ([20]) must be greater than zero unless all ui/s 
equal zero, which violates the condition that u>j ^ uik for 
j 7^ k. Thus, Equation (|2H|) does not hold, and d ^ 2. 
Moreover, since d = 2 is achieved by Algorithm 2, the 
theorem follows. ■ 



